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^ ■ Abstract 

h : 

This paper aims to present an elaborate view on the motivation 
and realization of the idea to extend Maxwell's electrodynamics to 
Extended Electrodynamics in a reasonable and appropriate way in or- 
der to make it possible to describe electromagnetic (3+l)-soliton-like 
objects in vacuum and in the presence of continuous media (external 
fields) [15]-[20], exchanging energy-momentum with the electromag- 
netic field. 
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1. Preliminary notes 

Classical Electrodynamics (CED) [1] is probably the most fascinating and 
complete part of the Classical Field Theory. Intuition, free thought, perspicu- 
ity and research skill of many years finally brought about the synthesis of 
experiment and theory, of physics and mathematics, which we have been 
calling for short the Maxwell equations for the duration of a century and a 
half. From the beginning of the second half of the 19th century till its end 
these equations turned from abstract theory into daily practice, as they are 
today. Their profound study during the first half of the 20th century brought 
forward a new theoretical concept in physics known as relativism. Brave and 
unprejudiced scientists enriched and widened the synthesis achieved through 
the Maxwell equations, and created a new synthesis known briefly as quan- 
tum electrodynamics. Every significant scientific breakthrough is based on 
two things: respect for the workers and their work and respect for the truth. 
"May everyone be respected as a personality, and nobody as an idol" one of 
the old workers used to say. We may paraphrase that saying: "may every 
scientific truth be respected, but no one be turned into dogma" . 

In our attempt to extend CED we tried to follow the values this creed 
teaches, as far as our humble abilities allow us to. Together with the analysis 
of the classical electrodynamics and the discrete conception for the structure 
of the electromagnetic field, the path followed brought us to the conclusion 
that a soliton-like solution of appropriate non-linear equations characterized 
by an intrinsic periodical process is the most adequate mathematical model of 
the basic structural unit of the field - the photon. The fact that neither the 
Maxwell equations nor the quantum electrodynamics offer the appropriate 
tools to find such solutions, unmistakably emphasizes the necessity to look for 
new equations. We decided not to follow the usual way for nonlinearization 
of CED [2], because, from our point of view, it does not comprise sufficiently 
fruitful and new ideas. Other attempts in this direction one may find in 
[3]-[12]. 

The leading physical ideas in our approach were the dual ("electro- 
magnetic") nature of the field on the one hand and the local Energy- 
Momentum Conservation Laws on the other. The realization that every 
such soliton-like solution determines in an invariant way its own scale fac- 
tor, as well as the suitable interpretation of the famous formula of Planck 
for the relation between the full energy E of the photon and the frequency 
v — 1/T of the beforementioned periodical process, which we prefer to write 
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down as h = E.T, helped us to formulate the rules for separating the realistic 
soliton-like models of the photon from among the rest. The resulting soliton- 
like solutions [19] possess all integral qualities of the photon, as described 
by quantum electrodynamics, but also a structure, organically tied to an in- 
trinsic periodical process, which in its turn generates an intrinsic mechanical 
angular momentum - spin (helicity). We consider this soliton- like oscillating 
non-linear wave much more clear and understandable than the ambiguous 
"particle-wave" duality. 

The dual 2-component nature of the field predetermined to a great extent 
the generalization of the equations in the case of an interaction with another 
continuous physical object, briefly called medium [20]. The proposed physical 
interpretation of the classical Frobenius equations for complete integrability 
[13], [14] of a system of non-linear partial differential equations as a criterion 
for the absence of dissipation, turned out to be relevant and was effectively 
used. The fruitfulness of the new non-linear equations is clearly shown by 
the family of solutions, giving (3+l)-dimensional interpretation of all (1+1)- 
dimensional 1-soliton and multisoliton solutions [20]. We note that making 
use of differential geometry proved extremely useful. 

This paper follows the main track of our efforts to build a clear and con- 
secutive picture of motivations and theoretical results. While most of the 
solutions to the new nonlinear equations we have found have been already 
published [15]- [20], a consecutive and well motivated reasoning, leading to 
the new equations is still missing. So, the stress in this paper is laid on the 
conceptual and generic framework. The purpose is to bring the reader to 
the conviction that the extension of CED, developed by us during the last 
10 years, is necessary and is physically and mathematically quite natural. 
We do not present solutions to equations. We set the problem to build a 
description of electromagnetic soliton-like objects, then we consider step by 
step physical arguments and try to find the corresponding step by step math- 
ematical adequacies. Our belief is that if the steps are in the right direction 
the positive results come inevitably. In its turn every new positive result 
(solution or relation), suggests new insights and invites to search for new 
more definite quantities and relations. This, step by step creative process, 
brought us to the today's motivational look on what we call Extended Electro- 
dynamics (EED). It worths from time to time to reconsider and re-estimate 
the importance and significance of any reason and idea been used, because 
this gives birth to new reasons and ideas and helps to sift out the basic and 
meaningful from the occasional and nonsignificant. 
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2. Physical conception for the electromagnetic field 
in Extended Electrodynamics 

As it is well known, the mathematical models of the real vacuum elec- 
tromagnetic fields in Classical Electrodynamics are "infinite", or if they are 
finite, they are strongly time-unstable. These models are not consistent with 
a number of experimental facts. A deeper analysis of CED, carried out in 
the first third of this century, resulted in the new conception for a discrete 
character of the field. This physical understanding of the field is the true 
foundation of EED and it shows clearly the principal differences between 
CED and EED. For the sake of clarity we shall formulate our point of view 
more explicitly. 

The electromagnetic field in vacuum is of discrete character and consists 
of single, not-interacting (or very weakly interacting) finite objects, called 
photons. All photons move uniformly as a whole by the same velocity 'c', 
carry finite energy 'E ', momentum 'p ' and intrinsic angular momentum. 
These features imply structure and internal periodic process of period T 
which may be different for the various photons. The quantity E.T', called 
"elementary action", has the same value for all photons and is numerically 
equal to the Planck constant 'h '. The invariance of 'c ' and 'h ' means nondis- 
tinguishability of the photons, considered as invariant objects. The integral 
value of the intrinsic angular momentum is equal to '±h '. For the topology 
of the 3-dimensional region, occupied by the photon at any moment, there 
are no experimental data, so it is desirable that the model-solutions to allow 
arbitrary initial data. 

We'd like to stress the following: EED considers photons as real finite 
objects, and not as convenient theoretical concepts, and it aims to build ade- 
quate mathematical models of these entities. So, the first important problem 
is to point out the algebraic character of the mathematical object describing 
a single photon. The corresponding generalization for a number of photons 
is then easily done. 

3. Choice of the modelling mathematical object 

According to the non-relativistic formulation of CED the electromagnetic 
field has two aspects: "electric" and "magnetic". These two aspects of the 
field are described by two 1-forms on TZ 3 and a parametric dependence on 
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time is possible: the electric field E and the magnetic field B. The following 
considerations will bring us to the conclusion that these two fields can be 
considered as two vector components of a new object, 1-form Q, taking values 
in a real 2- dimensional vector space, naturally identified with 1Z 2 . In fact, 
let's consider the question: do there exist constants (a, b, m, n), such that the 
linear combinations 

E' = aE + bB, B' = mE + nB 

give again a solution? The answer to this question is positive iff m — —b, 
n = a. The new solution will have energy density w' and momentum S' as 
follows: 

w> = i- (V) 2 + {B'f) = ^(a 2 + b 2 ) [e 2 + B 2 ) , 

S' = (a 2 + b 2 )^ExB. 

Obviously, the new and the old solutions will have the same energy and 
momentum if a 2 + b 2 = 1. 

This simple but important observation shows that besides the usual lin- 
earity, Maxwell's equations admit also "cross" -linearity, i.e. linear combi- 
nations of E and B of a definite kind define new solutions. Therefore, the 
difference between the electric and magnetic fields becomes non-essential. 
The important point is that with every solution (E, B) of Maxwell's equa- 
tions a 2-dimensional real vector space, spanned by the couple (E,B), is 
associated, and the linear transformations, which transform solutions into 
solutions, are given by matrices of the kind 

a b 
—b a 

If these matrices are unimodular, i.e. if a 2 +b 2 = 1, the energy and momentum 
do not change. It is well known that matrices of this kind do not change the 
canonical complex structure J in TZ 2 . (Recall that if the canonical basis of 1Z 2 
is denoted by (e±, e-i) then J is defined by J{e\) = e 2 , J{e-2) = — e±.) 

The above remarks suggest to consider E and B as two vector-components 
of an 7£ 2 -valued 1-form f2: 

Q = E <g> ei + B <g> e 2 . 
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So, the current j becomes 1-form J = j <g> e x with values in TZ 2 , and the 
charge density becomes a function Q = p <g> e 1 with values in TZ 2 . Maxwell's 
equations 

IdE Att 

_ rotB j, divB = 0, (1) 

c at c 

1 OB 

-— = -rotE, divE = Ann. (2) 
c dt ' w 

take the form 

1 <9Q 47r 

-^r = ^T-*dJ(fi), <^ = 4ttQ, (3) 

c at c 

where J(Q) = E <g> J(ei) + 5 ® J(e 2 ) = -E ® e 2 - 5 ® e x and d is the 
exteriour derivative. Note that according to the sense of the concept of 
current in CED and because of the lack of magnetic charges, it is necessary 
to exist a basis in TZ 2 , in which J and Q to have components only along e\. 
Nevertheless, this point of view shows that even at this non-relativistic level 
the separation of the EM-field to "electric" and "magnetic" is not adequate 
to the real situation. The mathematical object Q unifies and, at the same 
time, distinguishes the two sides of the field: there is a basis in TZ 2 , in which 
the electric and magnetic components are delimited, but in an arbitrary basis 
the two components mix (superimpose), so the difference between them is 
deleted. 

In the relativistic formulation of CED the difference between the electric 
and magnetic components of the field is already quite conditional, and from 
invariant-theoretical point of view there is no difference. However, the 2- 
component character of the field acquires a new meaning and manifests itself 
differently. In order to come to this we make the following considerations. 

As we mentioned above, some linear combinations of the electric and 
magnetic fields generate a new solution to Maxwell's equations. In particular, 
the transformation, defined by the matrix 

1 
-1 

defining a complex structure in TZ 2 , transforms a solution of the kind (E, 0) 
into a new solution of the kind (0, E) and a solution of the kind (0, B) into 
a solution of the kind {—B, 0), i.e. the electric component into magnetic and 
vice versa. This observation draws our attention to looking for a complex 
structure J, J 2 = —id in the bundle of 2-forms on the Minkowski space M, 
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such that if F presents the first component of the field, then J(F) to present 
the second component of the same field. Such complex structure really exists, 
in fact, it coincides with the restriction of the Hodge *-operator, defined by 
the pseudometric i] through the equation 



a A *(3 = —r](a, (3)y\r]\dx A dy A dz A = —r](a, (3)ol> , £ = ct, 

to the space of 2-forms: ** 2 — ~ ^a 2 (m)- So, the non-relativistic vector com- 
ponents (E,B) are replaced by the relativistic vector components (F,*F). 
The following considerations support also such a choice. 

The relativistic Maxwell's equations in vacuum dF = 0, d * F = are, 
obviously, invariant with respect to the interchange F — > *F. Moreover, if F 
is a solution, then an arbitrary linear combination aF + b * F is again a so- 
lution. More generally, if (F, *F) defines a solution, then the transformation 
(F, *F) — > (aF + b* F, mF + n* F) defines a new solution for an arbitrary 
matrix 

" a m 
b n 

Now, using the old notation VL for the new object f2 = F <g> e\ + *F <g> e 2 , 
Maxwell's equations are written down as dil = 0, or equivalently 5il = 0, 
where the coderivative operator 5 = *d* is just the (minus) divergence: 



-\7 U F". Clearly, an arbitrary linear transformation of the basis 



(ei,e2) keeps H as a solution. 

Recall now the energy-momentum tensor in CED, defined by 



-F aF a/3 5 u - F F V(J 



1 



in L 



(4) 



It is quite clearly seen, that F and *F participate in the same way in Q^, 
and the full energy-momentum densities of the field are obtained through 
summing up the energy-momentum densities, carried by F and *F. Since 
the two expressions F^F 1 "* and (*F) fia (*F) L ' a are not always equal, the 
distribution of energy-momentum between F and *F may change in time, 
i.e. energy-momentum may be transferred from F to *F, and vice versa. So 
we may interpret this phenomenon as a special kind of interaction between 
F and *F, responsible for some internal redistribution of the field energy. 
Now, in vacuum it seems naturally to expect, that the energy- momentum, 
carried from F to *F in a given 4- volume, is the same as that, carried from 
*F to F in the same volume. However, in presence of an active external 
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field, exchanging energy-momentum with f2, it is hardly reasonable to trust 
the same expectation just because of the specific structure the external field 
(medium) may have. So, the external field (further any such external field will 
be called also medium for short) may exchange energy-momentum preferably 
by F or *F, as well as it may support the internal redistribution of the 
field energy-momentum between F and *F, favouring F or *F. From the 
explicit form of the energy-momentum tensor it is seen that the field may 
participate in this exchange by means of any of the two terms F^ a F va and 
(*F) /ia (*F)' /a . Moreover, for the divergence of the energy-momentum tensor 
we easily obtain 



47T 



F lu/ (6F) v + (*F) IM/ (6*Fy 



(5) 



It is clearly seen that the quantities of energy-momentum, which any of the 
two components F and *F may exchange in a unit 4-volume are invariantly 
separated and given respectively by 

F^{5F)\ (*F) IM/ (8*F) V . 

But, in CED the exchange through *F is forbidden, the expression 
(*F)^ U (5*F) U is always, in all media, equal to zero. This comes from the un- 
conditional assumption, that the Faraday's induction law is universally true. 
Of course, we do not reject the existence of media, not allowing energy- 
momentum exchange through *F, but we do not share the opinion that all 
media behave in this same way. On the other hand, in case of vacuum, we 
can not delimit F from *F, these are two solutions of the same equation and 
it is all the same which one will be denoted by F (or *F), i.e. CED does 
not give an intrinsic criterion for a respective choice. Only in regions with 
non-zero free charges and currents, when dF = and SF = 4nj ^ 0, the 
choice can be made, but this presupposes (postulates) that the field is able 
to interact, i.e. to exchange energy-momentum, only with charged particles, 
i.e. through F. This postulate we can not assume ad hoc. 

Having in view these considerations we assume the following postulate in 
EED in order to specify the algebraic character of the modelling mathemat- 
ical object: 

In EED the electromagnetic field is described by a 2-form Q, defined on 
the Minkowski space-time and valued in a real 2-dimensional vector space W 



8 



and such, that there is a basis (ei,e 2 ) ofW in which Q takes the form 

n = F®d + *F®e 2 . (6) 

Since W is isomorphic to TZ 2 , further we shall write only TZ 2 and all 
relations obtained can be carried over to W by means of the corresponding 
isomorphism. In particular, every W will be considered as being provided 
with a complex structure J, so, the group of automorphisms of J is defined. 
Our purpose now is to prove that the set of 2-forms of the kind (6) is stable 
under the invariance group of J. 

First we note, that the equation aF + b * F = requires a = b = 0. In 
fact, if a ^ then F = - From aF + b*F = we get a*F-bF = and 
substituting F, we obtain (a 2 + b 2 )*F = 0, which is possible only if a = b = 
since *F ^ 0. In other words, F and *F are linearly independent. Let now 
(hi, k-i) be another basis of TZ 2 and consider the 2-form \l/ = G® k\ + hi- 
We express (k±, k 2 ) through (ei, e 2 ) and take in view what we want: 

G <E> ki + *G <S> k 2 = G (g> (aei + 6e 2 ) + *G <g> (me 1 + ne 2 ) = 

= (aG+m*G)<S>e 1 + (bG+n*G)<S>e 2 = (aG+m*G)®e 1 + *(aG+m*G)®e 2 . 

Consequently, bG + n*G = a*G — mG, i.e. (b + m)G + (n — a) * G = 0, 
which requires m = —b, n = a, i.e. the transformation matrix is 



a 


-b 


b 


a 



Besides, if Q\ and Q 2 are of the kind (6), it is easily shown that the linear 
combination XQi + is of the same kind (6). These results show that 
the 2-forms of the kind (6) form a stable with respect to the automorphisms 
of (1Z 2 ,J) subspace of the space A 2 (M, TZ 2 ). Moreover, if the component 
F is chosen, then the basis (ei,e 2 ) is uniquely determined. In other words, 
every fl of the kind (6) allows various representations of this kind, but every 
component F determines unique basis (ei(F) , e 2 (F)) , and so 2 subspaces 
{ei} and {e 2 } of TZ 2 . (Further the argument (F) of the corresponding basis 
is omitted). 

In order to separate the class of bases we are going to use, first we recall 
the product of 2 vector valued differential forms. If $ and W are respectively 
p and q forms on the same manifold N, taking values in the vector spaces 
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Wi and W 2 with corresponding bases (ei,...,e m ) and k n ), and </? : 

W^i x W 2 — > W3 is a bilinear map into the vector space W%, then a (p + 
form <£>(<&, \&) on N with values in W3 is defined by 

y? ($, *) = ^ A ¥ <g> y>(ei, fcj). 

In particular, if Wi = W 2 and = 1Z, and the bilinear map is scalar (inner) 
product g, we get 

5>* A 

Let now X and F be 2 arbitrary vector fields on the Minkowski space 
M, f2 be of the kind (6), Q^ v be the energy tensor in CED and g be the 
canonical inner product in 1Z 2 . Then the class of bases we shall use will be 
separated by the following equation 



Q^xny = l - * g (i(x)n, *i(Y)n) . 

We develop the right hand side of this equation and obtain 

±*g(i(X)Sl,*i(Y)ri) = 

i * g(i(X)F <g> ei + i{X) * F <g> e 2 , *i(Y)F ® e x + *i(Y) * F <g> e 2 ) = 
(i(X)F A *i(F)F)^( ei , ei ) + (i(X)F A «(F) * F)g(e u e 2 ) + 
+ (i(I)*FA«(y)f) 9 (e 2) e 1 ) + (i(X) * F A *i(F) * F)g(e 2 , e 2 ) 



(7) 



1 

= - * 
2 



— x^y 
2 



F°g{e x ,e x ) + (*F)^(*F)^(e 2 , e 2 )+ 
1 



= --X' i Y v 
2 



F^ + (*F)^(*F)' 



+ (F^(*Fr u + (*F)^F-)g(e 1 ,e 2 

In order this relation to hold it is necessary to have 

g(e u ei) = 1, g(e 2 , e 2 ) = 1, #(ei, e 2 ) = 0, 

i.e., we are going to use orthonormal bases. So, the stability group of the 
subspace of forms of the kind (6) is reduced to 5*0(2) or U(l). So, in 
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this approach, the group 5*0(2) appears in a pure algebraic way, while in 
the gauge interpretation of CED this group is associated with the equation 
dF = 0, i.e. with the traditional and not shared by us understanding, 
that the EM-field can not exchange energy-momentum with any medium 
through *F. 

4. Differential equations for the field 

We proceed to the main purpose, namely, to write down differential 
equations for our object Q, which was chosen to model the EM-field. We 
shall work in the orthonormal basis (ei,e2), where the field has the form 
(6). The two vectors of this basis define two mutually orthogonal subspaces 
{ei} and {e 2 }, such that the space 7Z 2 is a direct sum of these two sub- 
spaces: 1Z 2 = {e{\ © {e 2 }. So, we have the two projection operators 
7Ti : 1Z 2 — > {ei}, 7r 2 : 1Z 2 — > {e 2 }. These two projection operators extend to 
projections in the 7£ 2 -valued differential forms on M: 

7rifi = ir^n 1 <g) ki + Q 2 <S> k 2 ) = Q 1 <S> ^\k x + Q 2 <g> i\ x k 2 = 

= Q 1 ® 7Ti(aei + be 2 ) + fi 2 <S> ~K\{me\ + ne 2 ) = (aQ 1 + mVL 2 ) <g> e 1 . 
Similarly, 

7r 2 fi = (bQi + nfi 2 ) ® e 2 . 
In particular, if Q is of the form (6), then 

7Ti(F ® ei + *F <g> e 2 ) = F ® e ± , vr 2 (F ® d + *F ® e 2 ) = *F <g> e 2 . 

Let now our i?M-field Q propagates in a region, where some other contin- 
uous physical object (external field, medium) also propagates and exchanges 
energy-momentum with fl. We are going to define explicitly the local law 
this exchange obeys. 

First we note, that the external field is described by some mathematical 
object (s). From this mathematical object, following definite rules, reflecting 
the specific situation under consideration, a new mathematical object Ai is 
constructed and this new mathematical object participates directly in the 
exchange defining expression. The EM-field participates in this exchange 
defining expression directly through Q, and since Q takes values in TZ 2 , then 
Ai must also take values in 1Z 2 . 
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We make now two preliminary remarks. First, all operators, acting on 
the usual differential forms, are naturally extended to act on vector valued 
differential forms according to the rule D — > D x id. In particular, 

i i i i 

i i 

Second, in view of the importance of the expression (5) for the divergence 
of the CED energy-momentum tensor, we give its explicit deduction. Recall 
the following algebraic relations on the Minkowski space: 

5 P = (— l) p * -1 d* = *d*, 5* p = *d p for p = 2k+l, 5* p = — *d p for p = 2k. 

(8) 

If a is a 1-form and F is a 2-form, the following relations hold: 

* ( a A *F) = -cfF^dx" = * [(*F) A *(*«)] = \{*Fy v (*a) ^dx" . (9) 
In particular, 

*(F A *dF) = l-F^(dF)^ a dx a = *[6 * F A *(*F)] = (*F)^(5 * F) v dx» . 
Having in view these relations, we obtain consecutively: 
V„Q£ = V, [\F aP F^5; - F ucr F ua 

= ^F a ?V u F a/3 5; - {V v F, a )F™ - F^VvF™ = 
= l -F^[(dF) a ^ - VA - VpF m ] - {V V F^)F™ - F^V V F™ = 

= ^{dF) a ^ - F^V V F™ = -(*F)^V CT (*Fr - F^V a F™ = 

= (*F) lw (8*F) v + F lu/ (8Fy. 

Let now our field Q interact with some other field. This interaction, i.e. 
energy-momentum exchange, is performed along 3 "channels". The first 2 
channels are defined by the 2 (equal in rights) components F and *F of 
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f2. This exchange is real in the sense, that some part of the i?M-energy- 
momentum may be transformed into some other kind of energy-momentum 
and absorbed by the external field or dissipated. Since the two components 
F and *F are equal in rights it is naturally to expect that the corresponding 
2 terms, defining the exchange in a unit 4-volume, will depend on F and *F 
similarly. The above expression for V uQ 1 ^ gives the two 1-forms 

F^{5F) v d^, (*F)p,(6*F)"da? 

as natural candidates for this purpose. As for the third channel, it takes into 
account a possible influence of the external field on the intra-field exchange 
between F and *F, which occurs without absorbing of energy-momentum 
by the external field. The natural candidate, describing this exchange is, 
obviously, the expression 

F^(5 * F) W + (*F) fll/ (5F) l/ dx^. 

It is important to note, that these three channels are independent in the 
sense, that any of them may occur without taking care if the other two 
work or don't work. A natural model for such a situation is a 3-dimensional 
vector space K, where the three dimensions correspond to the three exchange 
channels. The non-linear exchange law requires some ^-valued non-linear 
map. Since our fields take values in TZ 2 this 3-dimensional space must be 
constructed from TZ 2 in a natural way. Having in view the bilinear character 
of V U Q^ it seems naturally to look for some bilinear construction with the 
properties desired. These remarks suggest to choose for K the symmetrized 
tensor product Sym(1Z 2 <g> 1Z 2 ) = 1Z 2 V 1Z 2 . So, from the point of view of the 
EM-Held, the energy-momentum exchange term should be written in the 
following way: 

*V(5tt,*tt). (10) 
In fact, in the corresponding basis (ei, e 2 ) we obtain 

* V (SQ, *Q) = * V (SF <g) e 1 + 5 * F <g> e 2 , *F <g> e x + * * F <g> e 2 ) = 

= *(5FA*F)(g)eiVei + *(5*FA*>KF)(g)e 2 Ve2 + *(-5FAF+5*FA*F)(g)eiVe 2 
= F^(5F) u dx» ® ei V ei + {*F^)(6 * Ffdx^ ® e 2 V e 2 + 
+ [F^(5 * F fdx» + (*F (U ,)6F"dx li ] <8> a V e 2 . 
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This expression determines how much of the EM-G.e\d energy-momentum 
may be carried irreversibly over to the external field (the first and the second 
terms) and how much may be redistributed between F and *F by virtue of 
the external field influence in a unit 4-volume. 

Now, this same quantity of energy-momentum has to be expressed by new 
terms, in which the external field "agents" should participate. Let's denote by 
$ the first agent, interacting with 7iiQ, and by ^ the second agent, interacting 
with 7T2^l. Since the corresponding two exchanges are independent, we may 
write the exchange term in the following way: 



According to the local energy-momentum conservation law these two quan- 
tities have to be equal, so we obtain 



This is the basic relation in EED. It contains the basic differential equa- 
tions for the _EM-field components and requires additional equations, spec- 
ifying the properties of the external field, i.e. the algebraic and differential 
properties of $ and \I>. The physical sense of this equation is quite clear: 
local balance of energy-momentum. 

The coordinate free written relationship (12) is equivalent to the following 
relations: ($ = a 1 <g> ei + a 2 <g> e 2 , ^ = a 3 <g> e 3 + a 4 <g> e 4 ) 



V ($,*7Tifi) + V(#,*7T 2 fi). 



(11) 



V (Stt, *n) = V($,*7Tifi) + V(^,*7T 2 fi). 



(12) 



SF A *F = a 1 A *F, 5*FA**F = a 4 A**F, 
5FA**F + 5*FA*F = a 3 A**F + a 2 A*F, 



(13) 




(14) 



Moving everything on the left, we get 

(SF - a 1 ) A *F = 0, (5 * F - a 4 ) A * * F = 0, 
(SF - a 3 ) A * * F + (5 * F - a 2 ) A *F = 0, 



or in components 



F^SF - a l f = 0, (*F)^(S * F - a 4 f = 0, 
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F^(5 * F - a 2 Y + (*F)^(5F - a 3 y = 0. 



Summing up the two equations 



F, U {5FY = F, u {a l y, (*F)^(6 * F) v = (*F)^(a A ) 



we obtain 



F^(6Fy + (*F)^(6*Fy = V 



V 



Ql = F^a 1 )" + (*F), u (a 4 y. 



This relation shows that the sum 



F^(a l y + (*F)^(a A ) 



is a divergence of a 2-tensor, which we denote by —P%- In this way we obtain 
the local conservation law 



Thus, we get the possibility to introduce the full energy-momentum tensor 



where is interpreted as interaction energy-momentum tensor. Clearly, P^ 
can not be determined uniquely in this way. 

So, according to (14), for the 22 functions F^ u , (a 1 )^ we have 12 equations, 
and these 12 equations are differential with respect to F^ v and algebraic 
with respect to Our purpose now is to try to write down differential 

equations for the components of the 4 currents a 1 . The leading idea in 
pursuing this goal will be to establish a correspondence between the physical 
concept of non- dissipation and the mathematical concept of integrability of 
Pfaff system. The suggestion to look for such a correspondence comes from 
the following considerations. 

Recall from the theory of the ordinary differential equations (or vector 
fields), that every solution of a system of ordinary differential equations 
(ODE) defines a local (with respect to the parameter on the trajectory) group 
of transformations, frequently called local flow. This means, in particular, 
that the motion along the trajectory is admissible in the two directions: we 
have a reversible phenomenon, which has the physical interpretation of lack of 
losses (energy- momentum losses are meant). Assuming this system of ODE 
describes fully the process of motion of a small piece of matter (particle), we 



v,(g; + p;) = o. 



(15) 



t; = q; + p;, 
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assume at the same time, that all energy-momentum exchanges between the 
particle and the outer field are taken into account, i.e. we have assumed that 
there is no dissipation. In other words, the physical assumption for the lack 
of dissipation is mathematically expressed by the existence of a solution - local 
flow, having definite group properties. The existence of such a local flow is 
guaranteed by the corresponding theorem for existence and uniqueness of a 
solution at given initial conditions. This correspondence between the mathe- 
matical fact integrability and the physical fact lack of dissipation in the simple 
case "motion of a particle" , we want to generalize in an appropriate way, 
having in view possible applications in more complicated physical systems, 
in particular, the physical situation we are going to describe: interaction of 
the field Q with some outer field, represented in the exchange process by the 
four 1-forms a 1 . This will allow to write down equations for a 1 in a direct 
way. Of course, in the real world there is always dissipation, and following 
this idea we are going to take into account its neglecting as conditions (i.e. 
equations) on the currents a\ As it is well known, the mathematicians have 
made serious steps towards study and formulation of criteria for integrability 
of partial differential equations, so it looks unreasonable to close eyes before 
the available and represented in appropriate form mathematical results. 

Remark. It is interesting, and may be suggesting, to note the following . 
In physics we have two universal things: dissipation and gravitation. We 
are going now to establish a correspondence between the physical notion of 
dissipation and the mathematical concept of non-integrability. As we know, 
the mathematical non-integrability is measured by the concept of curvature. 
General theory of Relativity describes gravitation by means of Riemannian 
curvature. The circle will close if we connect the universal property of any 
real physical process to dissipate energy-momentum with the only known so 
far universal interaction in nature, the gravitation. 

First we note that our base manifold, where all fields and operations 
are defined, is the simple 4-dimensional Minkowski space. According to our 
equations (12) the medium reacts to the field Q by means of the two 1Z 2 - 
valued 1-forms $ = a 1 ® e\ + a 2 <g> e 2 and \I/ = a 3 <8> t\ + a 4 <g> e 2 . So, we 
obtain four 7?.-valued 1-forms a 1 , a 2 , a 3 , a 4 . Because of the 4-dimensions of 
Minkowski space it is easily seen that only 1-dimensional and 2-dimensional 
Pfaff systems may be of interest from the Frobenius integrability point of 
view. All Pfaff systems of higher dimension are trivially integrable. 
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The integrability equations for 1-dimensional Pfaff systems are 

da* A a* = 0, i = 1,2,3,4. (16) 

Every of the 4 equations (16) is equivalent to 4 scalar nonlinear equations 
for the components of the corresponding 1-form. We note, that the solutions 
of (16), as well as the solutions of the general integrability equations for a 
p-dimensional Pfaff system, are determined up to a scalar multiplier, i.e. if 
a 1 define a solution, then fi.ct 1 (no summation over i), where f\ are smooth 
functions, define also a solution. 

In case of 2-dimensional Pfaff systems (a*, a- 7 ), defined by four 1-forms, 
their maximal number is 4.3=12. The Frobenius equations read 

da 1 Aq'Ao j = 0, i ^ j. (17) 

We have here 12 nonlinear equations for the all 16 components of a 1 , 
i = 1,2,3,4. Clearly, these equations (17) are substantial only if the cor- 
responding a 1 , the exteriour differential da 1 of which participates in (17), 
does not satisfy (16) or is not zero. 
Our assumption now reads: 

Every 2-dimensional Pfaff system, defined by the four 1-forms a 1 is 
completely integrable. 

As we mentioned above, physically this assumption means that we neglect 
the dissipation of energy-momentum. Note also, that 1-dimensional noninte- 
grable Pfaff systems are admissible, which physically means, that if there is 
some 1-dimensional nonintegrable Pfaff system among the a\ e.g. a 1 , then 
the corresponding dissipation of energy-momentum does not flow out of the 
physical system and it is utilized by the exchange processes, described by the 
rest currents a 2 , a 3 , a A . 

Finally we note, that (14) and (16) give in general 24 equations for the 22 
functions F^ v , (a) 1 , which seem to be enough to obtain the dynamics of the 
system at given initial condition. As for the Maxwell's theory, it corresponds 
to n 2 SQ = and ttiSQ = 47r(J/ ree + Jbound) <S> e±, i.e. a 2 = a 4 = and 
a 1 = a 3 = 4ir(J free + J bound )- 

In conclusion we note, that our notion of _EM-field requires a simultane- 
ous consideration of a number of soliton-like solutions and, in particular, a 
possible interaction (interference) among them. This problem does not fit to 
the goal we set in this paper, so it will be considered somewhere else. 
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